Numerical simulation of a multi-level atom interferometer 
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We present a comprehensive numerical simulation of an echo-type atom interferometer. The 
simulation confirms a new theoretical description of this interferometer that includes effects due to 
spontaneous emission and magnetic sub-levels. Both the simulation and the theoretical model agree 
with the results of experiments. These developments provide an improved understanding of several 
observable effects. The evolution of state populations due to stimulated emission and absorption 
during the standing wave interaction imparts a time-dependent phase on each atomic momentum 
state. This manifests itself as an asymmetry in the signal shape that depends on the strength of 
the interaction as well as spontaneous emission due to a non-zero population in the excited states. 
The degree of asymmetry is a measure of a non-zero relative phase between interfering momentum 
states. 



I. INTRODUCTION 

In recent years, atom interferometers (AIs) have be- 
come an invaluable tool for a variety of experiments re- 
lated to precision measurements and inertial sensing us- 
ing cold atoms [H4l3l]. 

The echo-type AI used in this work 0, Il3l - fl7| func- 
tions on the basis of phase modulation of the atomic 
wave function due to the interaction with standing wave 
(sw) pulses. Although this phase modulation is directly 
connected to the recoil energy, its functional form has a 
complicated dependence on a number of mechanisms — 
such as the dynamic population of magnetic sub-levels, 
phase shifts due to spontaneous and stimulated processes 
and the excitation of multiple momentum states. Some 
of these mechanisms have been studied in previous work 
EH i but many aspects of the AI are not fully un- 
derstood. In this work, we address all of these effects on 
the basis of a comprehensive numerical simulation and 
an improved analytical model. 

The simulation presented here successfully models data 
using measured experimental parameters as inputs and 
completes our understanding of a broad class of observ- 
able effects. One such effect is an asymmetry in the shape 
of signals produced by the AI, which is a manifestation 
of a non-zero relative phase difference between interfering 
momentum states. We show that the level of asymmetry 
is related to the amount of spontaneous emission occur- 
ring during sw excitations, as well as stimulated processes 
such as Rabi flopping. 

The analytical model of the AI accounts for effects due 
to both magnetic sub-levels and spontaneous emission. 
By including these effects, we avoid the need for a phc- 
nomenological model that was previously used in mea- 
surements of the atomic recoil frequency |16| . Using this 
theoretical treatment, we are also able to estimate mag- 
netic sub- level populations in the experiment. 

Although the analytical model provides an improved 
understanding of the response of the AI, the model is 
still limited to short pulses such that the motion of the 
atom during the interaction can be ignored (Raman-Nath 



regime). Additionally, the model is limited to detunings 
(A) large compared to the Rabi frequency (|A| » fl ) 
and the spontaneous emission rate (|A| 3> T). The an- 
alytical calculation also assumes that the excited state 
adiabatically follows the ground state, eliminating any 
dynamic exchange of amplitude or phase between states. 
The decay of the excited state due to spontaneous emis- 
sion is also accomplished in an approximate manner, 
since the ground state amplitude is not re-populated by 
the excited state. As a result of these limitations, the the- 
ory fails to model experimental data accurately in several 
regimes of interest. 

Since the simulation numerically solves the equations 
governing the system, none of the aforementioned limi- 
tations apply. The simulation accounts for the motion of 
atoms and the dynamic evolution of magnetic sub-levels 
populations during the interaction with sw pulses, which 
is not a feature of the analytical calculation. In general, 
this approach allows for a much broader class of phenom- 
ena to be studied. Additionally, a precision measurement 
of the recoil frequency using this technique would involve 
a detailed study of systematic effects. This justifies the 
need for an accurate and robust model for the signal. 

The rest of the paper is organized as follows. Section 
III1 gives a brief description of the experiment. Section 
IIIII reviews the main results pertaining to a theoretical 
calculation of the signal in two temperature regimes — 
one reminiscent of Bose-Einstein condensate (BEC) con- 
ditions involving one sw pulse, and the other similar to 
conditions in a magneto-optical trap (MOT), involving 
two sw pulses. A detailed description of the calculation 
of the one-pulse recoil signal is given in the Appendix. 
In £|IIIl we also derive an expression for the signal that 
accounts for magnetic sub-levels, which is used to model 
experimental data. In mVl we describe the theoretical 
background for the simulations. A discussion of the main 
results of this work follows in <JV] We present our con- 
clusions in ^VII and discuss future directions and appli- 
cations of this work. 
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FIG. 1. (Color online) Recoil diagram for a time-domain atom 
interferometer. Center of mass momentum states are shown 
as red dots. The 1 st sw pulse, applied at t = 0, diffracts 
the momentum states of the atom into multiples of the two- 
photon recoil momentum 2hk. The 2 nd sw pulse, applied at 
t = T, splits the momentum states further. Only the th order 
and the ±l st order diffractions from each sw pulse are drawn 
for simplicity. In the vicinity of the echo time, t = 2T, there 
is interference between all orders of momentum states. Our 
detection scheme is only sensitive to interferences between 
states that differ by 2hk. Two pairs of interfering momentum 
states are shown as solid black lines. 



II. DESCRIPTION OF EXPERIMENT 

The AI is used to measure the two-photon atomic re- 
coil frequency oj q = hq 2 /2M. Here, M is the mass of 
the atom and hq = 2hk is the momentum transferred 
to the atom from counter-propagating laser fields with 
wavelength A and wavenumber k = 2tt/X. A standing 
wave laser pulse with off-resonant traveling wave compo- 
nents interacts with a sample of laser cooled 85 Rb atoms 
(temperature T ~ 100 /iK) at times t — and t = T . 
During each pulse, atoms complete several two-photon 
transitions corresponding to the absorption of a photon 
from one traveling wave component and stimulated emis- 
sion into the oppositely directed traveling wave. This re- 
sults in the diffraction of atoms into a superposition of 
momentum states separated by hq, as shown in Fig.[TJ 

In the absence of spontaneous processes, and assum- 
ing both traveling wave components of the standing wave 
have the same polarization, the interaction returns the 
atoms to the same ground state magnetic sub-level. The 
atomic wave function develops a phase modulation on a 
time scale r q = 7r/w g ~ 32 /is, where r q is referred to 
as the recoil period. This phase modulation evolves into 
a spatial and temporal modulation in the atomic den- 
sity. Due to the finite velocity distribution of the atomic 
sample, the density grating dephases on a time scale de- 
fined by the coherence time t co h ~ 2/qa v ~ 1 /is. Here, 
a v ~ 10 cm/s is the width of the velocity distribution. 
Since t co h <C r q for the conditions of our experiment, an 
echo technique is used to cancel the effect of Doppler de- 
phasing and measure the temporal modulation induced 
on the density distribution due to atomic recoil. A sec- 
ond sw pulse, applied at i = T, causes interference of mo- 
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FIG. 2. (Color online) Pulse timing diagram for the exper- 
iment. After each sw pulse, a modulation in atom density 
forms and then decays in a time t co ^ ~ 2/qa v due to Doppler 
dephasing. At time t — 2T, the grating echo forms as a re- 
sult of the interference between different momentum states. 
A traveling wave read-out pulse coherently back-scatters light 
from the grating at this time. The intensity of this light is 
detected by a photo-multiplier tube. 



mentum states in the vicinity of the echo time, t = 2T, 
resulting in a rephasing of the density grating as shown 
in Fig. [21 The grating contrast is measured by applying a 
traveling wave read-out pulse (with the same wavelength 
as the sw pulses) and detecting the intensity of the co- 
herently back-scattered light from the atomic cloud. The 
density grating has various spatial harmonics in integer 
multiples of q due to the different orders of interfering 
momentum states. However, due to the nature of Bragg 
diffraction, this detection technique is only sensitive to 
density modulation that has a spatial periodicity of A/2 
(spatial harmonic q) — the smallest grating spacing that 
can Bragg scatter light of wavelength A. Only interfering 
momentum states that differ by Hq can produce such a 
modulation. The integrated intensity of back-scattered 
light (referred to as the echo intensity) is proportional to 
the contrast of the density grating. The echo intensity 
is a periodic function of the pulse separation, T, with 
period r q . 

The experiment utilizes 85 Rb atoms cooled in a 
magneto-optical trap contained within a glass vacuum 
system. The trap loads ~ 10 9 atoms in ~ 300 ms. Both 
the MOT laser beams and the magnetic field gradient 
are switched off prior to the AI experiment. The MOT 
B-field gradient is turned off in ~ 100 /Lts, after which 
the atoms are cooled in an optical molasses for several 
milliseconds. Three pairs of coils, one pair along each 
direction, are used to cancel residual magnetic fields and 
magnetic field gradients. These coils remain on continu- 
ously. Under these conditions, the magnetic field at the 
time of the echo experiment is canceled at the level of 
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1 mG over the volume of the trap. The AI excitation 
pulses are derived from two off-resonant, circularly po- 
larized, traveling wave beams. They are overlapped at 
the location of the trap to form a standing wave along 
the vertical direction. The AI beams have a Gaussian 
intensity profile and are collimated to a diameter of ~ 2 
cm. All excitation pulses are produced by using digital 
delay generators to trigger acousto-optic switches. The 
back-scattered light intensity from the grating is detected 
using a gated photo-multiplier tube. Measurements of u> q 
are accomplished by measuring the contrast of the grat- 
ing as a function of the pulse separation, T, on a suitably 
long time scale. A more detailed description of the ex- 
periment is given in Ref. [l5j . 



III. THEORY 

The theoretical expression for the recoil signal was 
originally derived in Ref. Q ■ The influence of magnetic 
sub-levels was subsequently included for studies of nano- 
structures in cold atoms [3]. The treatment in Ref. [l6| 
addressed the role of spontaneous emission on the recoil 
signal. In this work, we include both effects to develop 
a complete understanding of the properties of the recoil 
signal. Since theoretical derivations of this signal have 
been studied previously, we present the details in the ap- 
pendix. Here, we review the main results required to 
understand the recoil signal in the one-pulse regime (ap- 
propriate for BEC conditions) and the two-pulse regime 
(used for experiments under MOT conditions). These re- 
sults are used as a foundation for developing numerical 
techniques and for testing the accuracy of simulations. 



A. One-pulse recoil signal 

The one-pulse regime assumes the velocity distribution 
of the sample is infinitely narrow. The calculation is car- 
ried out in three stages. In the first stage, the Schrodinger 
equation is solved for the ground state amplitude of the 
atomic wave function, a g , using the two- level Hamilto- 
nian for a sw laser pulse of duration t% . The Hamiltonian 
assumes that the sw pulse is short (Raman-Nath regime) 
such that the motion of the atoms along the axis of the 
standing wave can be neglected during the interaction. 
This allows us to ignore the kinetic energy term in the 
Hamiltonian during the sw pulse. It also assumes that 
the pulse is far off-resonance (|A| 3> OqjT) such that 
the excited state is not significantly populated. Here, 
A = ui — luq is the detuning of the laser frequency, uj, 
from the atomic resonance frequency, uq, Oq = l^egEo/h 
is the Rabi frequency, fi eg is a dipole matrix element, 
Eq is the electric-field amplitude of each traveling wave 
component of the sw and 7 = T/2, where V is the spon- 
taneous emission rate. 

In the second stage of the calculation, after the pulse 
has turned off, the atom is allowed to evolve in free space 



for a time t, which results in a modification of the phase of 
the ground state amplitude. In an experiment involving 
a BEC, a traveling wave read-out pulse with wavelength 
A can be applied to the atomic sample and the inten- 
sity of the back-scattered electric field can be detected 
as the signal. The amplitude of the back-scattered elec- 
tric field is proportional to the A/2-periodic component 
of the atomic density modulation (the 2k = q Fourier 
harmonic) produced by the sw interaction. 

The final stage of the calculation requires a compu- 
tation of the probability density of the ground state, 
p g (r,t). The recoil signal is obtained by evaluating the 
square of the g-Fourier harmonic of this probability den- 
sity. 

The Hamiltonian in the field-interaction representation 
1191 for a two-level atom is 



n = n 



n(r) 







(1) 



where £l(r) = cos(fc • r) for a sw laser field. The en- 
ergy is defined to be zero for the ground state and — HA 
for the excited state. The — ihq term is a phenomcnolog- 
ical constant added to account for spontaneous emission 
during the interaction, which gives rise to amplitude de- 
cay of the excited state. This approach is valid in any 
open two-level system, but only approximately accounts 
for spontaneous emission in a closed system since the 
excited state population is not fed into the ground state 
(normalization is not preserved) . In addition, this Hamil- 
tonian does not account for the atomic recoil due to the 
jhoton emission. By using the density matrix approach 
2C I or resorting to Monte-Carlo wave function techniques 
2l| , a more complete model of spontaneous emission can 
be realized. 

As shown in the appendix, the amplitude of the back- 
scattered electric field as a function of the time, t, after 
the sw pulse is 

Ei(t) oc -uisin(y g t-0)[j o (/si) + J 2 («i)], (2) 



where u\ is the magnitude of the sw pulse area (Eq. lA4l) . 
uj q is the two-photon recoil frequency, 6 is a parameter 
associated with spontaneous emission given by 



= tan 



2AJ 



(3) 



J v (x) is the v th order Bessel function of the 1 st kind and 
Ki is given by Eq. lAT3l From Eq.|2]it is clear that the field 
amplitude is proportional to u\ in the small pulse area 
regime, and is shifted by a phase 9 due to spontaneous 
emission. This has interesting consequences for the signal 
shape and will be discussed in detail in flVB II 

For sufficiently cold atomic samples (~ 10 nK), t co h ^> 
T q and Doppler dephasing is negligible thereby allowing 
the temporal modulation in the contrast to be resolved 
after a single sw pulse [f| . The ground state density after 
the interaction with the sw pulse fEq. lAI0|) is shown in 
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Figs. 3(a) and |3(b)| for low and high pulse areas, respec- 
tively, in the absence of spontaneous emission (9 = 0). 
For small pulse areas, the density is sinusoidally modu- 
lated in space with a period A/2. The contrast of this 
spatial modulation oscillates in time with a period r q , 
as shown in Fig. 3(a) This density distribution corre- 
sponds to the momentum state \p = 0) interfering with 
the \p — ±hq) states. For higher pulse areas, the density 
distribution becomes more complicated (as shown in Fig. 
|3(b)[ ) due to the presence of additional spatial harmonics 
from the interference of higher order momentum states. 
The temporal modulation of the contrast also exhibits 
higher order harmonics but still primarily oscillates at 
the fundamental frequency 2u> q . 

The amplitude of the back-scattered field, E\ (t) (Eq. 
[2]) , derived from the ground state density p g (z, t) given by 
Eci lATOl is shown in Figs. |3(c)| and |3(d)| for the same pulse 
areas. The field changes with time due to the oscillations 
in the contrast of the density distribution. For small 
only the lowest order momentum states contribute to 
the signal, and Ei(t) oscillates sinusoidally at frequency 
uj q . In general, the states \nhq) and \(n+l)hq) inter- 
fere to produce a A/2-spatial modulation with a contrast 
that oscillates in time with a frequency 2(2n+ l)oJ q . For 
higher pulse areas, additional harmonics in the density 
distribution give rise to the shape of the back-scattered 
field shown in Fig. |3(d)| Equation [5] takes into account 
the interference of any two p-states differing by hq in mo- 
mentum, but no higher order interferences — for example, 
between \nhq) and \(n + 2)hq). 

The analytical expression for the one-pulse recoil sig- 
nal, denoted by s±, is simply the intensity of the back- 
scattered field: Si(i) cx |£i(£)| 2 . This quantity is propor- 
tional to the contrast of the density modulation. 



B. Two-pulse recoil signal 

The theoretical expression of the recoil signal in the 
two-pulse regime is a simple extension of that in the one- 
pulse regime. Typically the two-pulse experiment is car- 
ried out under MOT conditions where Doppler dephas- 
ing becomes important because the velocity distribution 
of the atoms can no longer be approximated by a delta 
function. 

As shown in the appendix, the standing wave interac- 
tion is equivalent to a phase grating in diffractive optics. 
The calculation for the two-pulse recoil signal involves ap- 
plying an additional phase grating e 1 ® 2 cos (i' r ) ( associated 
with the second sw pulse, to the atomic wave function a 
time t = T after the first pulse (shown in Eq. lA9l) . The 
atom evolves in free space after the second pulse until a 
time t { J2 o = (A + 1)T, called the A th order echo time, 
where the Doppler phases from all velocity classes cancel. 
Here, A = ri/Q, where ij is the momentum difference be- 
tween interfering p-states from the first pulse (in units of 
hq) and £ is the equivalent quantity for the second pulse. 
When averaged over the velocity distribution, the atomic 
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FIG. 3. Density distribution, p g (z,t), and corresponding 
back-scattered electric field amplitude, E\(t), for an atomic 
sample after the interaction with a weak sw pulse (a, c) and a 
strong sw pulse (b, d). In plots (a) and (b), z is the distance 
along the sw field. Dark portions correspond to low density 
while light portions correspond to high density. p g (z, t) shows 
spatial modulation that is periodic at integer multiples of A/2, 
and shows temporal modulation at integer multiples of twice 
the recoil frequency, 2ui q , as the pulse area increases. Equa- 
tion IA10I was used to produce plots (a) and (b) with first 
pulse area iti = 0.5 and 1.5 respectively. In plots (c) and (d), 
Eq. [5] was used with the same respective pulse areas. The 
zeroes in the contrast of the density modulation correspond 
to the zeroes in Exit). The effects of spontaneous emission 
were ignored for all plots. 



density shows modulation with a non-zero contrast only 
in the vicinity of these echo times. At all other times, the 
density is uniform since the modulation has been washed 
out by Doppler dephasing. In the general case, A is a 
rational number [lJ,LL9|], but in our case we consider only 
the first order echo which occurs at time: t^L = 2T, cor- 
responding to A = 1 . The density of the atomic sample 
after the two sw pulses and at the first order echo time 
is depicted in Fig. [5] 

The amplitude of the back-scattered electric field 
about t — 2T can be shown to be 



(Ea(At,T)) ocE^At} 



"(At/tcoh) 2 



x ysin" (u q (T + At) - 0) 

4 1 

Jo{k2) + - ^2(^2) + 2^4(^2) 



(4) 



where T is the separation between sw pulses, At = t—2T 
is the time relative to the echo time, 112 is the magnitude 
of the second pulse area and 

k 2 = 2u 2 \J s\u{uj q (T + At) — 9) s\n{u q (T + At) + 9) . 

(5) 

The (• • • ) brackets in Eq. @] indicate that the back- 
scattered field has been averaged over the initial velocity 
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distribution, which is assumed to be M axwellian with 
a characteristic width a v = ^j2ksT /M , which is also 
equal to the most probable speed. The coherence time, 
< co h = 2/qa v , is the time about the echo time for which 
the back-scattered field is non-zero and is determined 
solely by a v . 

From Eq.|4j which is equivalent to Eq. 25 in Ref. [l6j |. 
it is clear that for small U2 the basic signal dependence 
on pulse separation, T, is a sinusoidal oscillation at fre- 
quency Lu q that is phase shifted by an amount 8 (Eq.[3]). 
This phase shift is associated with spontaneous emission. 
As the second pulse area becomes large, the third factor 
in Eq.|4]has a more significant contribution to the signal. 
This is associated with the interference of higher order 
momentum states, as discussed in £1111 Al 

The two-pulse recoil signal is the integrated intensity 
of the back-scattered light in the vicinity of the echo 
time. The Ai-integrated intensity is approximately pro- 
portional to the square of the T-dependent part of Ei 
from Eq.gl 



s 2 (T)~^sin 4 ( Wg (T + Ai)-#) 
4 1 

^o(«2) + ^ -M^) + gJ4(K2; 



(6) 



As a function T, this expression is periodic with fun- 
damental frequency 2u> q (period r 9 = ir/ujq ~ 32 /is). 
The signal, §2, is proportional to the grating contrast 
measured in the experiment, which is a manifestation of 
phase modulation in the atomic wave function following 
the interaction with two sw pulses. The response of the 
AI depends on the level structure of the atom and the 
coupling between these levels and the driving field. We 
now explore these effects in detail. 



C. Recoil signal including magnetic sub-levels 

So far, we have reviewed the theoretical expressions for 
the one and two-pulse recoil signals, including the effects 
of spontaneous emission, for a two-level atom. In the ex- 
periment we use 85 Rb, which is a multi-level atom. If 
only the F = 3 — > F 1 = 4 transition is considered, there 
are 2F + 1 = 7 ground state magnetic sub-levels and 9 
excited state sub-levels. These energetically degenerate 
sub-levels have a significant role on the response of the 
AI. A previous treatment of the two-pulse signal [l4| ac- 
counted for multiple atomic levels, but assumed the pop- 
ulation was equally distributed and ignored effects due 
to spontaneous emission. Here, we extend the theoret- 
ical model of the recoil signal discussed in the previous 
two sub-sections to include multi-level atoms with an ar- 
bitrary distribution of initial sub-level populations. 

The coupling strength between states \g) = \n g J g m g ) 
and |e) = \n e J e m e ) is determined by the dipole matrix 
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FIG. 4. Comparison of two-pulse recoil curves predicted by 
the two-level theory (Eq. [6] — gray curve) and the theory in- 
cluding magnetic sub-levels (square of Eq. llObj — black curve) 
for various 2 nd pulse durations, T2 (a: 150 ns; b: 200 ns; c: 250 
ns; d: 300 ns). Here, the m- level populations were assumed 
to be equally distributed among the 7 levels of the J g = 3 



ground state. Pulse parameters: A 

r = r N . 



ior N ; n = 2.51* 



element 



where n a 



fi eg = -e {e 
= —e(n f 



c 9i 

JJ\r a 



r\g) 



g u gl^m g q L m c 



(7) 



n e are the principal quantum numbers, J g , J e 
are the total angular momenta and m g , m e are the mag- 
netic sub-levels of the ground and excited states, respec- 
tively. In our case, J g = F = 3 and J e = F' = 4. The 
unit vector e qL represents the polarization of the laser 
field. Linear and circular polarization states are denoted 
by qL (<Zl = for linear and qL = ±1 for er* polariza- 
tions, respectively). The factor (n e J e \\r\\n g J g ) in Eq.[7] 
is the reduced matrix element associated with the radial 
part of the wave functions the magnitude of which is 
unimportant for this treatment and will be absorbed into 
the Rabi frequency, Qq — fJ, eg E /h. The factor Cm g \ L m e 
is the Clebsch-Gordan coefficient, which describes how 
strongly two states are coupled by the photon and de- 
pends on the particular transition. Since we are only 
concerned with electric dipole transitions, it is non-zero 
only for states that obey the selection rules: J e — J g + 1 
and m e = m g + qL- 

From Eq. it is apparent that each degenerate m- 
level interacts with a sw pulse (of a given polarization) 
with a different coupling strength — which is proportional 
to the Rabi frequency for each transition: C m 9 \ L m c ^o- 
In the experiment, this differential coupling causes the 
population of the m-levels to become unbalanced after 
the interaction with the sw pulse (optical pumping) . The 
degree of the imbalance is determined by fio, A and the 
pulse durations. In the analytical treatment that follows, 
optical pumping is not taken into account. We assume 
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the population of each m-level remains constant during 
the sw pulses. However, the numerical simulation to be 
discussed in EHVI includes optical pumping effects. 

The area of a given sw pulse, denoted by index j ' = 1 
or 2, is given by 



,( m g) 



n 



3 



2A 



r 

2A 



1 -1/2 



fJg 1 Je 

Trig qL rn. 



(8) 



where Tj is the pulse duration. For the one-pulse and two- 
pulse signals, the back-scattered field amplitude from the 



state | J g Trig) is 



E 



C, 



Jg 1 Je 

nig q L m c 



x sm(uj q t - 9) [J (ki) + J2 («i)] 



(E ( r a \n *\(u^ } ) 2 (ci % i L t) 



sin 2 (uj q (T + At) -6) 

4 1 
^0(^2) + 2^2(^2) + -^4(^2) 



(9a) 



(9b) 



Here, it is understood that u\ in Eq. IA13I for k\ (for a 
two-level atom) has been replaced by u^ 1 ^ for the multi- 
level case. Similarly, u 2 in Eq.[5]for k 2 has been replaced 

by u^ 11 ^ ■ The extra factor of (Cm g \ L m e ) in Eqs. 1 
arises due to the coupling of states \ J g m g ) and \J e m e ) 
by the traveling wave read-out pulse (assuming that the 
scattered field has the same polarization as the read-out 
pulse) . The total scattered field from the atom is propor- 
tional to the sum of the fields scattered by each m-level 
weighted by the final population of that level, II j g mg : 

^{t) cx^n^™^ 7 "^) (10a) 

m g 

(E 2 (T))cxJ2^J e m g {E^\T)). (10b) 



The one and two-pulse recoil signals are proportional to 
the square of the total back-scattered field amplitudes 
given by Eqs. fTOl (sj oc l-E^I 2 , for index j = 1,2). The 
form of Eqs. [TOl allows for interference between scattered 
fields from each m-level. This additional interference 
from magnetic sub-levels strongly affects the shape of the 
recoil signal. Figure 0] shows a comparison of two-pulse 
recoil signals predicted by the two-level theory (Eq. [6]) 
and the theory including multiple sub-levels (square of 
Eq. UObp for various second pulse durations, T2. The 
two-level theory predicts extra zeroes in the signal shape 
that are not observed experimentally. For the same set 
of pulse parameters, the multi-level theory does not pre- 
dict these extra zeroes due to the interference of back- 
scattered light from each magnetic sub-level. The multi- 
level theory models experimental data much more effec- 
tively than all previous models, as we will show in ijV Al 



IV. DESCRIPTION OF SIMULATION 

The main goal of the simulation is to compute the two- 
pulse recoil signal (scattered field intensity as a function 
of pulse separation, T) under the conditions of the ex- 
periment. In practice, the signal in the experiment is 
obtained from a macroscopic sample (typically ~ 10 9 
atoms) at a temperature of T ~ 100 /iK. The most basic 
assumption of the simulation is that the signal from the 
AI can, in principle, be produced by a single atom. In- 
stead of velocity-averaging over many atoms with a well- 
defined momentum, as in the theoretical calculation of 
the two-pulse recoil signal, the simulation solves for the 
time-evolution of a single atom (represented by a wave- 
packet in momentum space) with a Maxwcll-Boltzmann 
probability distribution of velocities corresponding to a 
given temperature. 

The theoretical framework of the simulation uses the 
Schrodinger picture to account for a number of physical 
effects, such as the motion of the atom at all times during 
the evolution; the bandwidth of sw pulses; Doppler shifts 
of momentum states; momentum state excitation for ar- 
bitrary sw pulse duration, field strength and detuning; 
population and coherence transfer between ground and 
excited states; population dissipation due to spontaneous 
emission; and optical pumping of magnetic sub-levels. 

We solve for the time-dependent state amplitudes of 
the atomic wave function. The computational cost of this 
problem scales as N p multiplied by the number of time- 
steps used, where N p is the number of discrete momen- 
tum states. This is much less computationally expensive 
than the Heisenberg approach, where an N p x N p density 
matrix must be computed at every time-step. However, 
the Heisenberg approach allows for a much broader class 
of problems to be addressed, such as the atomic recoil 
due to spontaneous emission or AT-atom effects such 
as collective excitation and emission of radiation. Sim- 
ilarly, a Monte-Carlo wave function (MCWF) approach 
in the Schrodinger picture would allow for a more com- 
plete model of spontaneous emission [2l|, J2I] with addi- 
tional computational cost. In other work |15| . we carried 
out MCWF simulations of sw pulses to study the excita- 
tion of momentum states and their relative populations 
for comparison with an experiment. In the case of the 
two-pulse recoil experiment, we found that the random 
photon recoil due to spontaneous emission has a negli- 
gible effect on the signal shape compared to the excited 
state population damping associated with spontaneous 
emission. As a result, the photon recoil associated with 
spontaneous emission is not taken into account in the 
treatment presented here. 

We now review the theoretical framework of the sim- 
ulation. We restrict our model to atoms with two mani- 
folds of energetically degenerate magnetic sub-levels and 
total ground (excited) state angular momentum J g (J e — 
Jg + 1). We label the states of the atom as \J g m g ) 
and \J e m e ). The Schrodinger equation is solved numeri- 
cally in momentum space for wave functions in the field- 
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interaction representation under the influence of an off- 
resonant, time-dependent standing wave potential. The 
Hamiltonian for this system in the rotating wave and 
dipole approximations is (l9j 



P 2 h 



(11) 



x ( \p - hk) <p\ + \p + hk) (p | ) (S+ + 5- ). 



Here, P z is the momentum operator along the z-axis; 
A = to — uiq is the detuning of the laser field from res- 
onance; Sj g = \J g ) (Jg | (S Je = \ J e ) (J e | ) is the ground 
(excited) state projection operator; f2(t) = fiEo£(t)/H is 
the (time-dependent) Rabi frequency with envelope func- 
tion £(t); \p ± hk) (p | are raising and lowering operators 
for the momentum state \p) in units of the photon mo- 
mentum; and Sg L (S~ L ) is a raising (lowering) operator 
proportional to the atomic dipole operator for laser field 
polarization state qL. These operators are defined as fol- 
lows 



S £l I J 9 m a) 
S~ L | J e m e ) 
| J e m e ) 



= C, 
= C„ 



m, II ™ g +<iL \ J em g + q L ) (12a) 

m e \l m a -q L \Jg m e- 1l) (12b) 



= S QL \ J a m g) = 0- 



(12c) 



Clebsch-Gordan coefficients, Clh 'f,, describe the cou- 
pling of different ground and excited state m-levels. The 
pulse envelope, £(t), is based on a sum of exponential 
functions of the form l/(l + e~ x ), where x = (t — io)/T r isc 
for the rising edge of the pulse at time t = to, x = 
— (t — £i)/T r isc for the falling edge at t = t\ and r r i se is 
the rise time. Typically, the rise time is set equal to the 
experimentally measured value of 20 ns. The duration of 
the pulse is defined as r = t\ — to- 

The total momentum space wave function, $(p, t), can 
be written as a linear combination of the 2 J g + 1 ground 
states | J g m g ) 



and the 2J e + 1 excited states \J e m e ): 



m g 

+ S A/. TO. (p.*) 1 4 



(13) 



e ""el j 



where a j g mg and /3j e TOc are the corresponding state am- 
plitudes. 

The main goal of the simulation is to find the state 
amplitudes that simultaneously satisfy the Schrodinger 
equation and the set of rate equations that describe 
spontaneous emission. The solution to the Schrodinger 
equation is computed using a combination of two meth- 
ods. When the interaction potential is turned off (i.e. 
£(t) = 0) and the excited state population (Hj e (t)) is 
zero, an analytical solution of the Schrodinger equation in 
free space is used to compute <I>(p, t). At all other times, 
the solution is computed using a fourth order Runge- 
Kutta routine [23]. To describe spontaneous transitions 
from the excited to ground levels, the numerical solver 
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FIG. 5. (Color online) Evolution of ground state magnetic 
sub-level populations during a sw pulse for (a) a + and (b) 
linear polarization states. The gray grid lines in each plot 
indicate where the sw pulse turns on and off. For circular 
polarization there is optical pumping towards the m g = 3 
level, while for linear polarization the population is pumped 



towards the m g 
T risc = 20 ns, A 



= level. Pulse parameters: t% 

8.36r N ;fi = 5r N ;r = r N . 



300 ns; 



also simultaneously satisfies the following rate equations 
for the m-level populations: 



11 J a m„ 



= -rn 



J e m e 



V ( r h 1 



qL m g +q L 



rn 



J e m g +qL 



These equations have a solution 

n Jgmg (t>)=n Jgmg (t) + (i-e- r ( t '-v) 



(14a) 
(14b) 



(15a) 
(15b) 



1 Je 

J ra g q L m g +q L 



n 



7 C m g +q L 



it)- 



>IL 



Typically, T = Tn = 3.76 x 10' s , which is the natural 
radiative rate of the 5 2 P3/2 excited state for 85 Rb [24j . 

Figure O shows the evolution of m-level populations in 
the ground level during a sw pulse for (a) circular and 
(b) linear polarizations. Here, the initial populations are 
equally distributed among the 7 sub-levels. Evidence of 
optical pumping is clearly visible after only t\ ~ 300 ns, 
which is the typical duration of sw pulses in the experi- 
ment. 

The initial ground state population, nj s (0), is set to 
unity and the code ensures the normalization of the wave 
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function is preserved for all times thereafter. The popu- 
lation of a particular manifold at any time, Tlj(t), is the 
sum over all m- level populations, H/ m (t), in that man- 
ifold and is related to the state amplitude, aj m (p, t), 
through the following expressions: 

(*) = E n ^ 9 ™« (*) = !- n '. (*)> ( 16a ) 

nj 9 m 9 (t)= / |aj,m>,t)| 2 *- (16b) 



To account for the finite velocity distribution of the laser 
cooled sample, we use a Gaussian wave packet to describe 
the initial momentum distribution of each magnetic sub- 
level. The initial wave function can be written as 



$(p,0)=5>, /3m >,0)|J 9 



a j g mg (p, 0) = a Jg mg (0) 



,-(p-pof/2^ 



2tt(j v 



1/2 



(17a) 



(17b) 



Here, po is the average initial momentum, a p = ft/A is the 
standard deviation of the p-space probability distribu- 
tion and A = h/i^TrMksT) 1 / 2 is the thermal de Broglie 
wavelength. The coefficients aj m (0) are free param- 
eters that define the initial population of each ground 
state m-level. 

For the one-pulse recoil signal, the back-scattered field 
from each m-level is computed using the amplitude of the 
q- Fourier harmonic from each ground state amplitude, as 
given by 



Jg 1 J C 



x / a *j a m„ (p,t)a.J g m g (P ~ hq, t)dp. 



(18) 



This expression is analogous to Eq. [9a| in the analytical 
model. The total back-scattered field is a sum over the 
scattered fields from each m-level 



^(0 = 5>r' \ t ) 



(19) 



No m-level population factors appear in Eq.QlJl as they 
do in the analytical equivalent (Eq. |10a[) . because the 
populations are intrinsically built into the state ampli- 
tudes, a,j g mg , that appear in Eq.1181 The one-pulse recoil 
signal is then S\{t) = |23i(t)| 2 , where t is the time after 
the sw pulse. From Eq. [121 it is evident that the shape 
of the recoil signal depends on the distribution of initial 
ground state amplitudes, aj m _(0). 

In the two-pulse regime we use a temperature of ~ 50 
//K, which is 5000 times larger than that used in the one- 
pulse regime, but a factor of ~ 2 smaller than the typical 
temperature in the experiment. This is done to reduce 
the computation time in the two-pulse regime. We use 
a similar method to compute the intensity of the back- 
scattered field as in the two-pulse experiment. In the 



vicinity of the grating echo the intensity is \Ei(2T + t')\ 2 , 
where t' is the time relative to t = 2T and E\ is given 
by Eq. [THl Therefore, we calculate the two-pulse recoil 
signal by integrating this quantity over t' , as given by 



S2(T) 



J \E X {2T- 



t')\ 2 dt'. 



(20) 



This quantity is analogous to the expression for the two- 
pulse recoil signal for a two-level system (Eq.|BJ| or for an 
m-level system (square of Eq. llObl) . 

To conclude the description of the simulations, we 
discuss the discretization of momentum space. In the 
two-pulse regime, the width of the wave function is 
Op ~ 5hk. The discretization size is then typically set 
to Ap ~ (Tp/250 ~ ftfc/50 for good p-space resolution. 
However, due to the reduced temperature in the one- 
pulse regime, the width of the wave function is much 
smaller (o p ~ O.Olftfc — much less than the Hk transferred 
by absorption or emission along z). Hence, we can use 
Ap = Hk since this is the smallest momentum transfer 
allowed by the sw interaction. 

The number of sw pulses and the p-space discretization 
size has a large effect on the computational cost of the 
simulation in the two regimes. Since the two-pulse regime 
requires an additional sw pulse, a smaller Ap, and the 
evaluation of the signal is carried out as a function of T, it 
is much more computationally demanding than the one- 
pulse regime. Typically, the computation time is ~ 12 
hours for the two-pulse regime, and ~ 30 seconds for 
the one-pulse regime on a computing network using 2.2 
GHz processors (SHARCNET). However, the underlying 
physics in both regimes is essentially the same and many 
features of the two-pulse recoil signal can be studied using 
one-pulse simulations. 



V. RESULTS AND DISCUSSION 

Here we discuss the two main results of this work: the 
effects on the recoil signal due to magnetic sub-levels 
and the momentum-dependent phase imprinted on the 
ground state by spontaneous and stimulated processes. 
We also present, for the first time, a detailed description 
of the mechanisms that contribute to the shape of the 
signal. 



A. Effects due to magnetic sub-levels 

Figure [B] shows data from the two-pulse recoil experi- 
ment for various second sw pulse durations, r%. The echo 
intensity was recorded as a function of the pulse sepa- 
ration, T, for circularly polarized sw pulses {\qi\ = 1). 
Particular values of (70, 86 and 98 ns) were chosen to 
illustrate a range of signal shapes typically observed in 
the experiment. Two separate fits to the data are shown, 
one based on the two-level model (Eq.|BJ| and the other on 
the multi-level model (square of Eq. ll0b[> . The two- level 
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theory is insufficient to accurately model experimental 
data, whereas the multi-level theory fits all aspects of the 
signal for a large range of shapes. In particular, the multi- 
level theory successfully models the asymmetry and the 
broad valleys between zeroes of the signal that occur as 
the area of the second pulse is increased, as shown in 
Fig. 6(c) Fits using the multi-level theory show a factor 
of ~ 10 improvement in the x 2 /dof compared to that 
of the two-level theory, which corresponds to a factor of 
~ 3 improvement in the relative uncertainty of the recoil 
frequency. Thus, the multi-level model is ideally suited 
for precision measurements of u q using this technique. 

The development of this multi- level model — which con- 
tains only measurable parameters — represents a signifi- 
cant improvement in our understanding of the AI since 
previous efforts [l6j relied upon a phenomenological 
model to fit experimental data. This model was unable to 
explain the underlying physical mechanisms that govern 
the signal shape. Additionally, the present model is far 
more successful in fitting the full range of signal shapes 
that can be generated in the experiment compared to the 
phenomenological model. 

Another complication in modeling the AI is that mag- 
netic sub-levels and the spatial intensity profile of the 
excitation beams can be shown to produce similar effects 
on the signal shape. In this work, the beam diameter 
(~ 2 cm) was larger than the diameter of the atomic 
cloud (~ 1 cm). Therefore, we were able to demonstrate 
conclusively that magnetic sub-levels played the domi- 
nant role on the response of the AI. 

The data shown in Fig. [5] validates the multi-level 
model developed in ijlll CI A distribution of populations 
among several ground state magnetic sub-levels smears 
out any extra zeroes in the signal that would occur if 
the system were optically pumped into a single state, for 
example \ J g J g ). This smearing is due to interference be- 
tween the coherently scattered light from each m-level. 
Additionally, the valleys between the zeroes in the sig- 
nal are significantly broadened. These effects are most 
prominent when the pulse areas {u\ and u?) are large, 
since (1) higher order momentum states are contributing 
to the signal from each m-level (resulting in the double- 
peaked shape) and (2) population imbalance in the m- 
levels is maximized. 

It is also possible to estimate the final m-level popu- 
lations, ILj m , from the multi-level fits to the data in 
Fig. [SJ Although we found the statistical errors in the 
populations from the fits to be relatively large, the esti- 
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FIG. 6. (Color online) Data from the two-pulse recoil experi- 
ment for different 2 nd pulse durations, T2 (a: 70 ns; b: 86 ns; c: 
98 ns). Data is fit to the two-level expression (Eq.(6]) shown as 
the blue dashed line, resulting in a x 2 /dof ~9x 10 -3 , where 
dof is the number of degrees of freedom. The data is also fit 
to the multi-level expression (square of Eq . I10b|) shown as the 



red solid line, resulting in \ /dof ~8x 10~ . The \ ls com- 
puted assuming equally weighted data points (uncertainty for 
each point equal to unity). Final m-level populations, II j m , 
were estimated from the multi-level fits and are tabulated in 
Table [I] Pulse parameters: detuning A ~ 50 MHz; inten- 
sity I ~ 50 mW/cm 2 ; polarization state = 1; 1 st pulse 
duration n = 300 ns. 



T 2 [ns] 


n 3 -3 


n 3 - 2 


n 3 _i 


n 30 


n 3 i 


n 32 


n 33 


70 


0.00 


0.03 


0.32 


0.48 


0.15 


0.02 


0.00 


86 


0.00 


0.03 


0.38 


0.32 


0.21 


0.05 


0.01 


98 


0.00 


0.00 


0.50 


0.32 


0.15 


0.02 


0.01 



TABLE I. Final m-level population estimates from multi-level 
fits to experimental data for various 2 nd pulse durations, as 
shown in Fig.(5] The fits were performed assuming qL = 1. 



mates are considered accurate for two reasons. Firstly, a 
variation between the II j mg on the order of ~ 10% re- 
sults in a deviation in the minimum of the x 2 - Secondly, 
the distribution of populations is similar to those inferred 
from simulations. We attribute the large statistical errors 
to the presence of 12 free parameters in the fit function 
(based on the square of Eq. ll0a|) . namely u 2 , 8, At, u q , 
an amplitude factor A, and the 7 populations, ILj mg . 
From Table HI the general trend in the final populations 
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FIG. 7. Data from the two-pulse recoil experiment (a, c) 
and corresponding simulations (b, d) for two laser polariza- 
tion states (gray: ql = 0; black: |g_t| = 1) and two 2 
pulse durations (a, c: T2 = 150 ns; b, d: T2 = 250 ns). The 
recoil signals are qualitatively different for the two polariza- 
tions because of effects due to magnetic sub-levels. In the 
simulation, we used initial m-level populations If j g mg (0) = 



{0.05, 0.10, 0.20, 0.30, 0.20, 0.10, 0.05}, where m g 



-3 on the 



ment with the data for measured pulse parameters. Pulse pa- 
rameters: detuning A ~ 40 MHz; intensity / ~ 13 mW/cm 2 ; 
1 st pulse duration n = 500 ns; pulse rise time r r i sc = 20 ns. 



is to move towards the m g > levels as the pulse dura- 
tion increases. This is consistent with our expectations 
for c + excitation pulses. Magnetic sub-level populations 
can be confirmed using independent techniques [25j . 

Figure [7] shows data from the two-pulse recoil experi- 
ment obtained using circularly polarized beams, and sep- 
arately for linearly polarized beams. This data provides 
alternate confirmation of the role of magnetic sub-levels. 
The recoil curves for each polarization state are qualita- 
tively different for various second pulse durations. This 
figure also shows simulations that qualitatively agree 
with the data. 

The signal shapes shown in Fig.[7]for each polarization 
are consistent with our expectations based on two consid- 
erations. Firstly, the back-scattered field amplitude from 
each m-level is proportional to (C m 9 9 \ L m c ) IIj smg . For 
the same distribution of m-level populations at the time 
of the echo, the two polarizations of the read-out pulse 
induce different atomic responses and therefore produce 
distinct signal shapes for the two cases. Secondly, the 
atom-field couplings are larger toward the m g = J g state 
for the (?l = 1 polarization state (Cff | = 1) compared 
to the couplings for the m g = state (Cg \ f ~ 0.598). As 
a result, in the presence of spontaneous emission, circu- 
larly polarized sw pulses produces extreme state pump- 
ing. Similarly, for a linearly polarized sw pulse, the 

state 



atom-field couplings are largest for the m g 
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FIG. 8. (Color online) Data from the two-pulse recoil ex- 
periment for 2 nd pulse durations r-z = 70 ns (a) and 98 ns 
(b) overlayed with predictions of the two-level theoretical 
model (Eq. [6] — dashed blue curve) and a multi-level simu- 
lation (solid red curve). No fits were performed here. There 
is good agreement between the simulation and the data for 
measured pulse parameters. In contrast, the two-level the- 
ory gives a poor prediction of the data for the same pa- 
rameters. The distribution of initial m-level populations in 
the simulations was set to be similar to those shown Table 
11 n.j smg (0) = {0.05,0.10,0.25,0.35,0.20,0.03,0.02}. Pulse 
parameters: r\ = 300 ns; T r j so = 20 ns; A ~ 8.36 Tn; 
Q ~ 5.2 P N . 



however, the atom-field coupling is not as strong as it is 
for the m g = J g state with a er + -polarized beam. There- 
fore, for a given field strength, higher order momentum 
states will be populated for circularly polarized excita- 
tions. As a result of these two considerations, the recoil 
curves are qualitatively different for linear and circular 
polarizations. 

Figure [5] shows data overlayed with the prediction of 
the two-level model (Eq. [6J as well as output from nu- 
merical simulations for the conditions of the experiment. 
Only a temporal shift and amplitude scaling has been 
applied to the predictions of the analytical model and 
the simulation — no fits were performed in this plot. It is 
clear that the simulation successfully models the data for 
measured pulse parameters such as the detuning, pulse 
duration and intensity. The initial sub-level populations 
are not measured in the experiment and were assumed 
to be similar to those measured in Table U (distribution 
peaked near m g = 0). In contrast, the two-level model 
gives a poor prediction of the signal shape when the mea- 



sured pulse parameters are used as inputs. These results 
suggest that the simulation can be used to study a va- 
riety of properties associated with this interferometer, 
including possible systematics that would affect a pre- 
cision measurement of the recoil frequency. We discuss 
other applications of the simulation in WII 



B. Phase evolution during standing wave 
interaction 

We now describe the phase evolution of momentum 
states subject to standing wave excitation. We separate 
the discussion into three parts: phase dynamics due to 
spontaneous emission; phase dynamics due to stimulated 
emission and absorption; and a measurement of the com- 
bined effects including magnetic sub-levels. These effects 
are expected to be important for matter wave interference 
when the excited state contributes significantly to the 
relative phase between the different momentum states of 
the ground level. We describe these effects in the context 
of a single sw interaction using a theoretical treatment 
and compare the theoretical predictions to the results of 
simulations. Finally, we observe trends associated with 
these effects using the echo AI and verify these observa- 
tions using two-pulse simulations. 




t (us) 

FIG. 9. Comparison of recoil curves (back-scattered field 
intensity as a function of time) predicted by the one-pulse 
theory (Eq. IA16[) with and without spontaneous emission. 
The black curve corresponds to a spontaneous emission-free 
system (T = — s- 6 = 0), while the gray curve corre- 
sponds to a system where spontaneous emission is present 
(r = T N 9 = -0.133 rad). The solid grid line shows 
the first zero in the signal after t = in the absence of 
spontaneous emission, while the dashed grid line shows the 
zeroes shift by a temporal amount St = 8/ui q in the pres- 
ence of spontaneous emission. Pulse parameters: A = 7.5Tn; 
Qo — 1.5 Tn; Ti = 250 ns; pulse area ui ~ 1.43. 



terference by truncating the infinite sum in Eq. IA11I to 
only n = —2, . . . , 2 so that: 



1. Spontaneous emission 

In this section, we summarize our theoretical under- 
standing of the influence of spontaneous emission on the 
shape of the recoil signal. The expression for the one- 
pulse recoil signal (Eq. lA16|) contains the essential physics 
relevant for this discussion. The consequences of sponta- 
neous emission on the two-pulse recoil signal (Eq.|5|) stem 
from the same effects as those in the one-pulse signal. 

There are two main features that are observed as a con- 
sequence of including spontaneous emission in the theory. 
Firstly, within a single period, maxima adjacent to the 
signal zeroes have different amplitudes. We refer to this 
as the asymmetry in the recoil signal. Secondly, there is 
a temporal shift of the zeroes in the signal toward earlier 
times relative to the zeroes without spontaneous emis- 
sion. Both of these features are best demonstrated when 
the pulse area is large — causing the double-peak struc- 
ture within each period of the recoil signal. This is il- 
lustrated in Fig. [9] on the basis of Eq. lAl6l and has been 
observed experimentally (see Figs. |H1 [7] and Both 
of these features are present only when the spontaneous 
emission parameter, 9, is non-zero (see Eq.[3]). 

For the moment we consider a system without spon- 
taneous emission (8 — and the complex pulse area, 
0i = u\e %B , is real). The signal develops multiple peaks 
within each period for large pulse areas (compare Fig. 
3(c) to Fig. |3(d)| . This structure is due to the inter- 
ference of higher order momentum states which become 
more populated as u\ increases. We demonstrate this in- 



zSi(<)«-J 2 (0i)Ji(0 1 )e- j3 ^ t 
- Ji(0i)Jo(0i)e~^ 4 + J o (0i)Ji(0i)e^ 4 
+ Ji{e 1 )J 2 (Q* 1 )e l3 ^ t + J 2 (0i)J 3 (0t)e l5 ^ 4 . 



(21) 



This expression contains the first three harmonics of the 
back-scattered electric field that oscillate with frequen- 
cies ±uig, ±3uj q and 5uj q , respectively. The sum of differ- 
ent harmonics gives rise to constructive and destructive 
interference in the field amplitude. The result is a dimin- 
ished amplitude between the zero crossings of the total 
back-scattered field — producing the double-peaked struc- 
ture. This effect is shown in Fig. 10(a) using the terms 
corresponding to n = 0, 1, and 2 in Eq. [21] In the ab- 
sence of any phase shifts of the individual harmonics, the 
shape of the recoil signal (si(t) = \Ei(t)\ 2 ) is symmetric 
(black curve in Fig.|9|). 

We now consider the case where spontaneous emission 
is present in the system (0^0). The spontaneous emis- 
sion parameter, 9, causes a phase shift of each harmonic 
comprising Ei(t) in Eq.O This is demonstrated by ex- 
panding the Bessel functions in Eq.[2l]in a Taylor series 
and keeping only the leading terms: 



E x {t) « - ( m + J sm(uj q t -0) + ^ sm(u q t + 6) 



+ ^ sin(w 9 t - 36) - ^ sin(3w 9 < - 6) 
8 8 

+ -i sin(3cV + 0) + -± sin(3uV - 30). 
o4 9d 



(22) 



12 



i 



0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 



■ (2hk\ihk) 
(4hk\6hk) \ 






All interferences 



10 20 



30 40 
T (us) 



50 60 



1 



0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 







\ / 1 


All interferences ■ 



10 20 30 40 50 



60 



T (us) 



FIG. 10. (Color online) Real part of the harmonics comprising 
the back-scattered electric field, Ei(t) (Ea JAll)) . for a system 
without spontaneous emission (a: 9 — 0) and with sponta- 
neous emission (b: 9 = —0.133 rad). Each graph shows har- 
monics associated with the interference of momentum states 
\nhq) and \(n + l)hq) for n = 0,1 and 2 (red, green and 
blue curves, respectively). These harmonics correspond to 
the n = 0, 1 and 2 terms in Eq. lAllI The total back-scattered 
field, shown as the black curve, is the sum of all harmonics 
due to the interference between p-states differing by hq. In 
(b), the total scattered field is shifted by a temporal amount 
St = 9/uj q as shown by the dashed grid line. Note that the 
imaginary part of E\ (t) sums to zero in both (a) and (b) when 
all terms in Eq. lAllI are included. The pulse parameters are 
the same as in Fig.(9] 



Each term in Eq.[2U has been phase shifted by an integer 
multiple of 9, but no two terms are shifted by the same 
temporal amount (colored curves in Fig. [10(b)] ). When 
all values of n are included, the different temporal shifts 
sum to a net shift of St — 9/oj q (shown by the dashed 
grid lines in Figs. [5] and |10(b)[ ). 

Figure |10(b)| illustrates the harmonic components of 
the back-scattered field amplitude (E\ (t) from Eq.[2]) as- 
sociated with the interference (nhq |(n + l)hq) for n = 0, 
1 and 2 using 9 = —0.133 rad. When summed over all 
harmonics, the field exhibits more constructive interfer- 
ence on one side of the waveform than the other (black 
curve in Fig. |10(b)[ ) giving rise to an asymmetry in the 
peak amplitude within each recoil period and a tempo- 
ral shift of the zeroes. The corresponding field intensity, 
si(t), is shown in Fig. [S] as the gray curve. When 6 = 0, 
there is no temporal shift, thus no asymmetry, and we 
recover the black curve in Fig. [5] 

In the experiment, we always observe a negative asym- 



metry, which means that the amplitude of the first peak 
of the recoil signal is higher than the second peak be- 
tween consecutive zeroes, as shown in Figs. [S] and [51 A 
positive asymmetry corresponds to the second peak am- 
plitude being larger than the first. A negative asymmetry 
implies that the zeroes of the signal are shifted towards 
earlier times (9, St < 0), which is a consequence of using 
a positive detuning. These observations are consistent 
with theoretical predictions. 

The phases of each harmonic in Eq. [22] represent the 
difference in phase between two interfering momentum 
states. This phase difference is a direct result of the de- 
cay of the excited state into the ground state. When 
the harmonics are summed over, the relative phases are 
effectively averaged and the resulting phase shift of the 
scattered field amplitude represents a weighted average 
phase difference between all interfering momentum states 
that differ by hq. To first order, this phase is approxi- 
mately equal to 9. Spontaneous emission is evidently one 
mechanism by which the relative phases of two interfering 
momentum states can be non-zero. In the next section, 
we will show that stimulated emission and absorption 
produce similar effects. 



2. Stimulated emission and absorption 

In this section, we illustrate that stimulated emission 
and absorption between the ground and excited levels 
leads to a time-dependent asymmetry in the recoil sig- 
nal. To isolate this effect, we generate simulations of 
the one-pulse recoil signal in the absence of spontaneous 
emission (i.e. F = 0). Under these conditions, the signal 
exhibits either positive or negative asymmetry depend- 
ing on the pulse duration, t\, as shown in Fig. 1111 This 
originates from the dynamic evolution of the atomic state 
amplitudes during the sw pulse (Rabi oscillations) due to 
stimulated emission and absorption. In contrast, the an- 
alytical theory for the one-pulse recoil signal (Eq. lA16l) . 
which ignores the dynamic exchange of population be- 
tween the ground and excited states (see Eq. lA2ap . pre- 
dicts that there should be no asymmetry or phase shift 
when there is no spontaneous emission. Results from the 
simulations indicate that, by changing t\, the asymme- 
try oscillates between positive and negative at the same 
frequency as the generalized Rabi frequency. 

The oscillations in the asymmetry shown in Fig. [11] 
originate from the time-dependence of the relative phase 
between interfering momentum states. As the interac- 
tion time with the sw pulse increases, there is a trans- 
fer of momentum from the ground state to the excited 
state and back to the ground state again. The interme- 
diate process of going through the excited state imprints 
a momentum-dependent phase on the ground state wave 
function. As a result, the phase difference between two 
adjacent p-states changes as a function of the interaction 
time, T\. 

The amplitude of the phase oscillation is largely de- 
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FIG. 11. (Color online) Simulations of one-pulse recoil sig- 
nal, si(t), in the absence of spontaneous emission for pulse 
durations varying within one Rabi oscillation (n ~ 100 — 116 
ns). In part (a) the recoil signal begins roughly symmetric 
at ri = 100 ns. As ri is increased the signal exhibits posi- 
tive asymmetry and shifts toward later times. Part (b) shows 
the corresponding Rabi oscillation in the ground state popu- 
lation for each n in (a). In part (c) the signal shifts toward 
earlier times and shows negative asymmetry as ri increases 
from 110 to 116 ns. Part (d) shows the corresponding Rabi 
oscillations for (c). Here, the generalized Rabi frequency is 
5.4 Tn ~ 32 MHz. Pulse parameters: Trjse ~ 1 
fio = 2 Tn; r = 0. 



ns; A = 5 Tn ; 



pendent on the bandwidth of the sw pulse, which has 
contributions from the pulse duration, n, and the rise 
time, T r i sc . A large bandwidth pulse produces a large am- 
plitude phase oscillation. As T\ increases, the bandwidth 
of the pulse decreases — giving rise to a damping of the 
oscillation amplitude. This result can also be inferred 
from the time-dependence of the excited state popula- 
tion, since the relative phase between momentum states 
is a measure of its contribution to the recoil signal. 



3. Combined phase effects 

In the experiment, both stimulated and spontaneous 
processes contribute to phase shifts of the recoil signal. 
In the following section, we discuss the combination of 
these two mechanisms and measure their affect on the 
signal. This constitutes one of the most important results 
in this paper. 

We first describe an analytical calculation that predicts 
a dynamic phase for each momentum state that depends 
on both the parameter 6 (Eq. [3]), which is a measure 
of the number of spontaneous emission events, and the 
pulse area u\, which is a measure of the number of stimu- 
lated processes. This calculation is valid for the one-pulse 
recoil signal. 

Since the interferometer is sensitive to the relative 
phase of interfering momentum states, this p-dependent 
phase provides a direct connection between matter wave 



interference and transitions due to both spontaneous and 
stimulated processes. However, the calculation of the rel- 
ative phase is subject to the limitations of the theoreti- 
cal framework developed in £11111 (weak, far off-resonance 
pulses and short interaction times) . This necessitates the 
use of simulations to develop more accurate predictions. 

When spontaneous emission is included in the theory, 
the excited state imprints a momentum-dependent phase, 
4>{p), on the ground state wave function, which can be 
obtained from 



cj}(p) = arg[a g (p)], 



(23) 



where a g (p) is the amplitude of state \p) in the ground 
level. The phase of a discrete momentum state \p = nhq) 
can be obtained using Eq. lA7l for a g (p): 



4 1 )=arg[H)"J n (u 1 e ie )] i 



(24) 



where the superscript (1) indicates that this calculation 
is valid for a single sw pulse and the subscript n labels 
the momentum state. Expanding the Bessel function in 
a series in powers of u±e , one can show that 
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To first order: 4>n ~ \n\6 — rm /2, which is independent 
of the pulse area, Ui, and therefore ignores stimulated 
processes. This result indicates that the phase difference 
between interfering momentum states is 



= 0« - <j>. 



(i) 



+ 7r/2 for n > 1, 
+ 7r/2 for n< 1. 



(26) 



For higher orders, however, stimulated processes con- 
tribute to the phase of each momentum state and there 
is a distinct dependence on the pulse area. The AI is sen- 
sitive to the average relative phase between p-states that 
differ by hq. For an arbitrary ground state amplitude 



j(p), this quantity can be expressed as 



(Ac/)) = arg 



(27) 



where the (• • • ) notation denotes an average over mo- 
mentum. For the specific case of the one-pulse recoil 
experiment, the average relative phase is 



(Acj) {1) ) = arg 



Jn(uie lB ) J n -x(uxe ie ) 



, (28) 



which has a distinct dependence on the sw pulse duration, 

Tl. 
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FIG. 12. (Color online) Average relative phase between mo- 
mentum states differing by hq as a function of sw pulse du- 
ration, n. An offset of tt/2 has been added to the vertical 
axis. The dashed gray grid line shows the value of 6 ~ —0.06 
rad from Eq.[3] The dashed red line is a calculation based on 
Eq. [28] with u\ given by Eq, IA4l The dotted blue line shows 
the prediction of Eq. [27] computed from a simulation in the 
absence of spontaneous emission. The solid black line is the 
analogous prediction of a simulation including spontaneous 
emission. The steep rise of the two simulations at small n 
represents the first Rabi cycle of phase transfer from p — 
to neighboring states. These states initially have zero ampli- 
tude and their phases are not well defined. Pulse parameters: 
T risc ~ 1 ns; A = 8.4 E N ; flo = 4.5 Tn- 



Figure IT21 illustrates the dependence of (A^ 1 )) on the 
sw pulse duration, t\. Predictions are shown for both 
theory (dashed red line — Eq. l2"51) and simulations in the 
presence and absence of spontaneous emission (dotted 
blue and solid black lines, respectively — Eq. l27l) . 

For small t\ , the theoretical expression for the average 
relative phase (Eq .125)1 is approximately equal to 6. This 
result can also be inferred from Eq. l26[ which is valid only 
for small u\. According to theory, the contribution from 
spontaneous emission is a constant phase offset equal to 
9 (Eq. [3]). As the interaction time increases, the theory 
predicts an increase in |(A0W)|, as well as a slow os- 
cillation. The frequency of this oscillation is empirically 
determined to be ~ 251q/A ~ 2.7 MHz, which is due to 
stimulated emission and absorption. 

Results of simulations in the absence of spontaneous 
emission (dotted blue line) show an additional high fre- 
quency oscillation in (Ac^ 1 )) as a function of T\. These 
oscillations are caused by the transfer of phase between 
the ground and excited state due to stimulated processes, 
and have a frequency approximately equal to the gen- 
eralized Rabi frequency a/Hq + A 2 ~ 56 MHz. These 
high frequency phase oscillations are not taken into ac- 
count in Eq. [28j as a result of the far off-resonance as- 
sumption. When spontaneous emission is included in the 
simulations (solid black line), the oscillations are highly 
damped. As the pulse area increases, there is a signifi- 
cant departure of the theory from the simulations. We 
attribute this departure to the absence of dynamic pop- 
ulation transfer between the ground and excited states 
and the weak-field assumption made in the theory. 
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FIG. 13. (Color online) Variation of average relative phase, 
(A^ 2 '}, as a function of 2 nd pulse length. The relative phase 
measured from fits to experimental data (♦) are overlayed 
with the results of a two- level simulation (A, where J g — 0) 
and a multi-level simulation (•, where J g = 3). A subset of 
the multi-level simulations overlayed with data is shown in 
Fig- [HI (A^ 2 ') is measured from the data and the multi- 
level simulation by fitting to the square of Eq. llObl and using 
the best fit value of 9. Similarly, it is measured for the two- 
level simulation by fitting to Eq. [6] A subset of the data 
from which {A<^)' 2 ') is extracted is shown in Fig. [6] with the 
corresponding fits shown as solid red curves. The error bars 
for each set of points are equal to the statistical uncertainties 
in 6 from the corresponding fits. The pulse parameters are the 
same as those in Fig. [6] except for the two- level simulations 
where the rise time was reduced to r r i sc = 10 ns to illustrate 
the presence of oscillations in (Ac// 2 - 1 } as a function of T2. 



We now review results of measurements of the average 
relative phase, (Acj)^), from the two-pulse recoil experi- 
ment and corresponding simulations. Figure [T3"l shows the 
variation in (A</>( 2 )) as a function of t%. This variation 
corresponds to the change in the asymmetry of the sig- 
nal shapes shown in Fig. [5] We measure (A0( 2 )) from the 
data by fitting to the square of Eq. UObl (which includes 
effects due to m-levels) and using the best fit value of 9. 
From the data, (Aip^) seems to decrease and level out 
near —35 mrad after ti ~ 90 ns. This can be explained by 
the excited state population approaching steady state — 
giving rise to a saturation of momentum state excitation. 

Figure [13] also shows results of two-level simulations 
(J g = 0) implemented with the same pulse parameters 
as the data (except for the rise time which was set to 
Trisc = 10 ns). Here, (A</>( 2 )) is measured by fitting each 
recoil curve generated by the simulation to Eq.[6] (which 
incorporates only two levels). This is equivalent to the 
method used to measure the average relative phase from 
data. The trend of the two-level simulation is qualita- 
tively similar to the black curve shown in Fig. 1121 In 
comparison to the data we find an additional superim- 
posed oscillation and a constant phase offset. The fre- 
quency of this phase osc illation is the same as the general- 
ized Rabi frequency {\Jti^ + A 2 ~ 59 MHz in this case). 
This is consistent with effects due to stimulated emission 
and absorption, as discussed in WB 21 The oscillations 
are quickly damped as the excited state population ap- 
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proaches steady state. The rate at which this occurs is 
mainly determined by the bandwidth of the pulse (which 
has contributions from both t 2 and r r i se ) and the spon- 
taneous emission rate. In these two-level simulations, we 
have reduced the rise time of the sw pulses from 20 ns 
(as measured in the experiment) to 10 ns in order to 
demonstrate the presence of these oscillations. 

The red circles shown in Fig. [13] represent simula- 
tions including magnetic sub-levels (J g = 3). These re- 
sults show no phase oscillations as a function of inter- 
action time, which is consistent with the data. The ab- 
sence of phase oscillations in the data and the multi-level 
simulation is attributed to the relatively long rise time 
(irise ~ 20 ns) of the excitation pulses, as discussed in 
the previous sub-section. There is also additional damp- 
ing of phase oscillations due to the interference between 
back-scattered fields from each ground state magnetic 
sub-level. 

Although the trend in (Acj)^} from both the two- level 
and the multi-level simulations is similar to the data, 
there is a small phase offset between the two. This can 
be explained in the following manner. Firstly, the offset 
is equivalent to a difference in the level of asymmetry 
between the data and the simulations. This is not sur- 
prising since our estimates of various input parameters, 
such as the Rabi frequency and m-level populations, are 
not directly measured in the experiment. Secondly, our 
ability to measure (Acj)^) from the data is limited by the 
signal-to-noise. Since (Acf)^} is estimated by fitting to 
the data, and each point in the fit is equally weighted, this 
results in an over-estimation of the phase value and an 
under-estimation in the uncertainty of that value. This is 
especially true at small t%. This is because the quantity 
being measured is highly shape-dependent and, at small 
T2, where there is very little asymmetry in the data, a 
small amount of noise can cause a large error in the mea- 
surement. 

The time evolution of the relative phase predicted by 
the two-level and the multi-level simulations, however, 
are very similar. This suggests that (A</>( 2 )) is largely 
unaffected by the presence of magnetic sub-levels, but de- 
pends mostly on the sw pulse parameters and the magni- 
tude of the spontaneous emission rate — which, aside from 
the pulse rise time, were the same in both simulations. 
This is consistent with the theory presented above. 

In a previous paper (l6| . we investigated the influence 
of spontaneous emission on the two-pulse recoil signal us- 
ing both theory and experiments. We used a phenomeno- 
logical model for the two-pulse recoil signal that was suc- 
cessful in predicting the asymmetry in the signal shape 
and the temporal offset. Since the spontaneous emission 
parameter, is related to the spontaneous emission rate, 
r, it appeared that there was a variation in T with the 
pulse duration and a measurement of T appeared to be 
possible by fit ting experimental data using this model 
(Eq. 26 in Ref. |l6[). However, the average relative phase 
of the signal, (Atfi^), depends on both 9 and the pulse 
area. The phenomenological model did not take into 



account this pulse area dependence. As a result, erro- 
neous values of 9, and therefore T, were measured when 
T2 changed. Thus, the conclusion that V depended on T2 
was not justified. The dependence of the average relative 
phase on spontaneous and stimulated processes, as well 
as effects due to magnetic sub-levels, was not fully un- 
derstood until recently. With this interpretation, we can 
explain the variation of 9 with ti observed in Ref. [16]. 

VI. CONCLUSION 

We have developed a comprehensive numerical simula- 
tion of a time-domain echo interferometer that is success- 
ful in predicting a variety of physical effects associated 
with matter wave interference. The simulation accurately 
models experimental results for measured input parame- 
ters. Based on these results, we have determined that the 
interference between back-scattered electric fields from 
each ground state magnetic sub-level significantly affects 
the signal shape, which is consistent with experimental 
observations. Additionally, a more complete understand- 
ing of the phase modulation of the atomic wave function 
induced by standing wave pulses has been developed on 
the basis of the simulation. 

The simulation motivated the development of an im- 
proved analytical model of the signal that accounts for 
both spontaneous emission and ground state magnetic 
sub-levels. Using the analytical model, we were able to 
improve the relative error in measurements of the recoil 
frequency by a factor of ~ 3. We were also able to esti- 
mate magnetic sub-level populations at the time of the 
read-out pulse. 

Through theoretical analysis and simulations, we have 
clarified the role of spontaneous emission on the recoil 
signal. We have also explained how the relative phase 
between interfering momentum states affects the signal 
on the basis of stimulated emission and absorption, which 
is most important when the excited state is significantly 
populated. These effects are not unique to the AI con- 
sidered in this paper, but apply more generally to in- 
terferometers that involve momentum state interference 
induced by laser excitation, such as frequency domain 
AIs [26[. 

The robustness of the simulations suggest that they 
can be used to study the AI under a variety of interesting 
conditions, such as the near- resonance (| A| ~ T), high in- 
tensity (f2 ~ |A|) or long-pulse (r, >• |A| _1 ) regimes. 
An important application of this work also includes the 
study of systematic effects on a precision measurement 
of the atomic recoil frequency using the three-pulse tech- 
nique outlined in Refs. [Ill, Eg]- 
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APPENDIX 

We review the calculation of the recoil signal in the 
one-pulse regime. The evolution of the atomic state am- 
plitudes is governed by the Schrodinger equation 




(Al) 



where T-L is the Hamiltonian given by Eq.[TJ We assume 
that the rate of change of the excited state amplitude 
is small compared to the detuning (|d e | <C |A|), which 
results in the following relations: 



it) 



(A 2 



7 2)l/2 "9' 

'(A 2 + 7 2 )V2 a 



(A2a) 
(A2b) 



Here, 9 is a phase that characterizes the amount of spon- 
taneous emission occurring during the pulse and is given 
by Eq. [3] Equation IA2bl is integrated to obtain the 
ground state amplitude for a pulse of duration t\ 



a g (r) = a g (r,0)e 



-zOi cos q r 



(A3) 



where q = k\ — &2 = 2fc is the difference between the k- 
vectors comprising the sw and ®\ = u\e is the complex 
pulse area with magnitude u\ given by 
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(A4) 



We assume the initial state of the atom to be a plane 
wave, a g (r,0) = V^ 1 ' 2 e lp0 ' r ' ,i , momentum po = 

Mvo = hko and interaction volume V. 

In Eq. lA3l we have ignored a — i@i term in the exponent 
because it is independent of r and is therefore unimpor- 
tant for interference. It also causes a g to decay exponen- 
tially during the pulse (since 0i is a complex quantity), 
which does not preserve normalization. Since we assume 
the system is closed (no atoms are being lost) this corre- 
sponds to an unwanted physical process and is a direct 
result of using a Hamiltonian that is non-Hermitian (Eq. 

From Eq. IA3[ it is clear that the effect of interacting 
with the sw pulse is to spatially modulate the phase of the 



ground state amplitude. Equivalently, the standing wave 
field is a phase grating that diffracts the atomic wave 
function into a superposition of momentum states. Using 
the Jacobi- Anger expansion, we decompose the ground 
state amplitude into its harmonics: 



a g (r) = a g {r,Q) V (-»)» J^Q^e 1 



(A5) 



where J n (x) is the n th order Bessel function of the 1 st 
kind. Here, the index n indicates the order of each har- 
monic, which corresponds to a specific momentum state: 
e mq r ^ |p = r j,^g), xhe sw pulse only populates the 
ground state with even multiples of the photon momen- 
tum, hk. In contrast, the excited state is populated with 
only odd multiples of hk due to an extra factor of fl(r) 
in the excited state amplitude, a e , which can be written 
as: 



a e (r) 



2(A 2 + r 2 )V2 



a g (r,0) 



^(-i) n J„(6i) e^ 2 ™- 1 )^ + e i(2n+l)fe-r. 



(A6) 



The population of each momentum state, | J„(0i)| 2 , is 
completely determined by the magnitude of the pulse 
area, u\. However, the phase of each momentum state is 
determined by both the magnitude and the phase of the 
pulse area, u\ and 9, respectively. This has an important 
consequence for the interference of two momentum states 
(see CVBH) . 

After the sw pulse turns off, the atom evolves in free 
space. The population of the excited state is scaled by 
a factor f2 /|A| <§; 1 compared to the ground state, and 
is therefore ignored. The Hamiltonian in free space con- 
tains only the kinetic energy term, P 2 /2M, where P is 
the momentum operator. Therefore, the calculation of 
the ground state amplitude at time t is much more trans- 
parent by transforming into momentum space: 



a g(P) 



(2irh) 



3/2 



V 



The solution to the Schrodinger equation in this regime 
dictates that the phase of the ground state amplitude is 
modulated differently for each momentum state: 



,(p,i) = a ff (p)exp 



i p 
h2M 



(A8) 



By transforming back to position space, the ground state 
amplitude a time t after the sw pulse is: 



1 



a g (r,t) = —7=e 



i(po-t — eo*)/R 



V 



Y^{-i) n Jn{ < di)e mq ' r e~ mq ' Vat e 



(A9) 
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where eo = Po/2M is the initial kinetic energy of the 
atom, v = p /M is the initial velocity and to q = Hq 2 /2M 
is the two-photon recoil frequency. 

We are interested in how the ground state atomic den- 
sity distribution, p g (r,t) = \a g (r,t)\ 2 , evolves in time: 



■71,7) 



(A10) 



x Q-imr giW-fot e irj(2n+ri)uj q t 



This quantity contains a sum over all orders of interfer- 
ence between momentum states and is shown in Figs. 3(a) 
and 3(b) The index 77 — equal to the difference between 
momentum states in multiples of Hq — denotes the order 
of the interference. This form of p g helps identify the 
phase due to spontaneous emission and its subsequent 
role in the one-pulse recoil signal. The density is mod- 
ulated in time at integer multiples of u q . This multiple 
depends on the order of the momentum state, n, and the 
order of the interference, ij. The Doppler phase rjq ■ Vot 
results in a dephasing of the density modulation when 
the density is averaged over the velocity distribution of 
the atomic sample. However, for the purposes of this dis- 
cussion, we have assumed that the velocity distribution 
is infinitely narrow such that Doppler dephasing is negli- 
gible. As a result, we have ignored the Doppler phase in 
all expressions that follow. 

The density modulation can be detected in an exper- 
iment by measuring coherent Bragg scattering from a 
traveling wave read-out pulse. However, this technique 
is sensitive only to the g-Fourier harmonic of the density 
distribution when the read-out pulse has the same wave- 
length, A, as the sw pulse. This is because Bragg scat- 
tering occurs for structures with a spatial periodicity of 
to A/2, for m = 1,2,.... Higher-order spatial harmonics 
of the density grating (i.e those periodic at A/4, A/6, . . .) 
do not scatter the light. Thus, the back-scattered electric 
field amplitude is proportional to the q- Fourier harmonic 
of p g (r,t), which corresponds to the amplitude of the 



e lvq r term in Eq. lAlQI with 77 set to unity: 

fi(t)a J2 *Jn(ei)Jn+im)e l{2n+1) ^ t . (All) 



A simpler form of the back-scattered electric field am- 
plitude is realized by using the Bessel function summa- 
tion theorem [H G3, Sz| 



^j„(e 1 )j n+ ^(ej 



*\ i(2n+rj)iAi q t 



sin(tV -0) \ v/2 
sin(uj q t + 6) 



(A12) 



where 

K\ = 2ui\Jsin(u} q t + 8) sm(uj q t — 6), (A13) 
resulting in the following expression for E\ (t) 



Exit) cx -Ji(rti) 



sin(oj g t — 9) 
sm(oj q t + 0) 



1/2 



(A14) 



However, this form contains a singularity at t = — (6 + 
TO7r)/w 9 , which is not physical. The singularity can be 
removed by using the recurrence relation: 2vJ u {x) — 
x\J v —i(x) + J u+ i{x)\ , such that 

Ei(t) cx -uisin(LJ q t-6)[j (Ki) + J 2 («i)]. (A15) 

The analytical expression for the one-pulse recoil signal, 
denoted by si, is the intensity of the back-scattered field, 
which is proportional to \E\(t)\ 2 : 

h(t) oz u 2 sin 2 (Lo q t - e)^^) + J 2 (ki)] 2 . (A16) 

This expression for the recoil signal is valid for an atomic 
sample with a narrow velocity distribution (BEC condi- 
tions) after interacting with a single sw pulse. 
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